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Abstract 

A new algorithm for computing the accurate values of quasinormal fre- 
quencies of extremal Reissner-Nordstrom black holes is presented. The nu- 
merically computed values are consistent with the values earlier obtained 
by Leaver and those obtained through the WKB method. Our results are 
more precise than other results known to date. We also find a curious fact 
that the resonant frequencies of gravitational waves with multi-pole index 
I coincide with those of electromagnetic waves with multi-pole index I — 1 
in the extremal limit. 



1 Introduction 

The determination of the quasinormal frequencies of black holes has for several 
years been of interest and is related to experiments which aim at detecting 
gravitational waves from supernovae or coalescences of binary neutron stars, 
which are thought to eventually form a black hole. In the late stage of the process 
of black hole formation, a certain mode of the gravitational wave dominates the 
emission. This is called the quasinormal mode of a black hole. 

*TIT/HEP-311/COSMO-61 
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The equations governing oscillations of black holes are derived by Regge 
and Wheeler Q, and Zerilli|^] for the non-rotating uncharged case, Zerilli|| and 
MoncriefQ for the non-rotating charged case and TeukolskyQ for the rotat- 
ing uncharged case. The resonant frequencies of the Schwarzschild black holes 
were first computed by Chandrasekhar and Detweiler^jj, who treated it as a 
boundary value problem of the second order ordinary differential equation of 
Regge- Wheeler. Integration of the Regge- Wheeler equation for the quasinormal 
boundary condition is numerically unstable. To avoid the instability, LeaverJjJ 
presents the numerically stable continued fraction method for Schwarzschild and 
Kerr black holes. Here, a quasinormal mode function corresponds to a mini- 
mal solution of the recurrence relations which are satisfied by the coefficients 
of a series expansion of a wave function. The minimal solution of the three- 
term recurrence relation is obtained by the corresponding continued fraction^. 
The continued fraction method can give the values of frequencies with high 
numerical precision because it uses no approximation, as is the case in the semi- 
analytic WKB method developed in other works & [II], fll ^2|. Leaver ]l3[ also 
generalizes the continued fraction method to calculate accurate values of the 
quasinormal frequencies of Reissner-Nordstrom black holes, though Kokkotas 
and Schutz[^l| raised the question of the applicability of the continued fraction 
method to charged black holes. 

We are indeed motivated by his paper fll3|]. There is another interesting 
question, that is, whether a continued fraction method can be used for the 
maximally charged black holes not discussed in his paper. In the extremal limit, 
the wave equation has an irregular singular point at the horizon, which makes the 
series expansion of a solution there invalid, and thus a continued fraction method 
does not seem to work. In this paper we show that a continued fraction method 
is applicable even to the extremal black holes if we expand a solution about a 
suitable ordinary point. In section II, we derive an eigenvalue equation with 
continued fractions after expanding a wave function about an ordinary point 
to determine the quasinormal frequencies of the extremal Reissner-Nordstrom 
black hole. In section III, we summarize the numerical results and find an 
interesting coincidence between the resonant frequencies of gravitational waves 
with multi-pole index I and those of electromagnetic waves with multi-pole index 
Z — 1 in the limit of maximal charge. Section IV is devoted to conclusions and 
discussion. 



2 Eigenvalue Equation for Quasinormal Modes 

Leaver obtained the recurrence relations of Schwarzschild, Kerr[0, and Reissner- 
Nordstrom black holes [jl3| by expanding solutions about the event horizon and 
then solved the characteristic equations with a continued fraction to get these 
resonant frequencies accurately. However his method is not valid when a black 
hole is extremal. The wave equation has an irregular singular point at the 
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horizon about which the expansion of the solution is not available because a 
radius of convergence for a series expansion about an irregular singular point 
generally vanishes. Therefore we have to expand the solution about a suitable 
ordinary point in the extremal black holes. 

We start from the Zerilli-Moncrief equation for perturbations on Reissner- 
Nordstrom backgrounds of mass M and charge Q written as 
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and s = 0,1,2 are for a massless scalar field, an odd parity electromagnetic 
perturbation, and an odd parity gravitational perturbation, respectively. The 
new parameter p s is introduced here in order to deal with the case of the scalar 
field as well. Note that the definitions of q% and are exchanged with each other 
from the book of Chandrasekhar|T^] , though the most of other notations follow 
this book. The scaling of t and r in this paper is such that c = G = M = 1. 
We also use p = -iw as a frequency variable for convenience. 

Now we follow Leaver's procedure of getting a recurrence relation which 
plays a central role in a continued fraction method. The tortoise coordinate in 
the extremal limit becomes 

r, =r + log(r-l) 2 -^-, (11) 



which makes a difference in the basic equations between extremal and non- 
extremal black holes. The Zerilli-Moncrief equation (ffl) is therefore rewritten in 
terms of r-coordinate as follows, 
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This equation has irregular singular points at r = 1 and r = oo as we mentioned 
in the previous section, which makes the scries expansion about the horizon 
invalid because it is hopeless that the series expansion is convergent. Hence it 
is necessary to improve Leaver's treatment of obtaining the quasinormal modes 
of the extremal black holes. 

The boundary conditions of a quasinormal wave function are given by 

Z s ~ e~ pr * as r* -> oo, (13) 
Z s ~ e pr * as r* — > -oo, (14) 

which mean there only exists purely outgoing wave at the infinity and purely 
ingoing wave at the black hole horizon. These are also transformed to 

Z s ~ r- 2p e- pr as r -> oo, (15) 
Z s ~ e pr (r- l) 2 P e - p ^ as r -► 1, (16) 

which determine a pre-factor of a quasinormal wave function in the following 
form 

oc 

Z s = e pr e- p ^ (r - ifP r -*P e -Mr-i) a n u n , (17) 

n=0 

where u is specified by a point where we expand a solution. We use a series 
expansion about an ordinary point, which has a finite radius of convergence at 
least up to the nearest singular point. Thus, we use u = (r — 2)/r instead of 
u = (r — l)/r which is used in Ref. [|l3| . In our choice of u the black hole horizon 
and infinity correspond to u = — 1, 1, respectively. The distances from the point 
u = to two singular points u = —1,1 are equal, which can make us possible to 
examine the boundary conditions at both sides simultaneously. It is clear that 
the quasinormal mode boundary conditions are satisfied if and only if Ea„u™ 
converges at both boundaries, horizon and infinity. This is equivalent to the 
condition that both Ea„ and £(— ) n a n are finite, that is, both Sa2n and £<Z2 n +i 
are convergent simultaneously. To examine the convergence of each summation, 
we first obtain the recurrence relation of coefficients a n of the expansion (|l7|), 
and then attempt to obtain the recurrence relation only for b n — a2 n +i and that 
of c„ = am ■ The quasinormal mode boundary conditions are satisfied when the 
odd and even sequence are minimal solutions of the corresponding recurrence 
relations, respectively. 

After substituting the above series expansion (p~7]) into Eq.([ll) we obtain a 



five-term recurrence relation in the following, 

"102+7100 = 0, (18) 

"203 + 7 2 ai + 5 2 a = 0, (19) 

a n a n+ i + (3 n a n + 7„a„_i + <5„a„_ 2 + e n a n -3 = 0, (n > 3), (20) 
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where the recurrence coefficients are given by 

a n = n 2 +n, (21) 

Pn = 0, (22) 

7„ = -4A - 2n 2 - p s + 2q s + n{2 - 24p) + 12 p - 64p 2 , (23) 

S n = 2(p s -q s ), (24) 

e„ = n 2 -p s - 12p+ 16p 2 + n(-3 + 8p). (25) 



For the scalar field, we can find that (3 n — S n — in the above relation, which 
means an odd sequence and even sequence are completely decoupled with each 
other, because the Zerilli-Moncrief equation in terms of u is symmetric under 
the transformation u — > —u. Therefore it has only symmetric or anti-symmetric 
solutions as quasinormal mode functions. This is the same as the case of a 
quantum mechanical system with a reflection symmetry, where a wave function 
of energy eigenstate must be either symmetric or anti-symmetric. This makes 
the calculation of quasinormal frequencies of the scalar held easy. Each sequence 
is characterized by the three-term recurrence relation ( |20| ) starting from the 
two-term relation (|l|) or |[|). We can obtain two independent characteristic 
equations below similar to one for the quasinormal frequencies of Schwarzschild 
black holes, 

02 _ _7i_ 
a at 

03 _ _72 

These two equations are for the modes of even-numbered and those of odd- 
numbered overtone, respectively. 

For gravitational or electromagnetic quasinormal modes, we have to develop 
a little more complicated way because the presence of <5 n -term makes odd and 
even sequences dependent on each other, though it seems that they are not much 
correlated for large n due to the fact that the coefficient 6 n becomes relatively 
small as n increases compared with other terms. We first need to eliminate 
the odd or even sequence to examine the convergence of the summation of 
each sequence. This procedure is presented in Appendix A. For b n a five-term 



recurrence relation below is obtained, 

aife 2 +/3i6i+7i^o = 0, (28) 

a 2 b 3 + $ 2 b 2 +72&1 +S 2 b = 0, (29) 

UnK+l + Pnbn + InK-l + t>nK-2 + = 0, (n>3), (30) 

and for c n we get a following relation, 

<5 2 c 3 + f3 2 c 2 + 72C1 -I- S 2 c = 0, (31) 

+ p n c n + 7„c„_i + <5„c„_ 2 + £nC„-3 = 0,(n> 3), (32) 
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where each recurrence coefficient is also given in Appendix A. 

As mentioned previously, the quasinormal mode boundary conditions are sat- 
isfied if and only if both sequences are minimal solutions of recurrence relations 
simultaneously. In order to use a continued fraction method, both five-term 
recurrence relations are transformed to three-term recurrence relations by suc- 
cessive Gaussian elimination steps which are exhibited in Appendix B. Then, 
the ratios of the first two terms of both sequences for minimal solutions are 
determined by the corresponding continued fractions in the following, 



a 2 ex p»- p»- p»- 

03 = h = 

01 h 0»- p»- 



(33) 
(34) 



where the quantities with double-prime are recurrence coefficients of trans- 
formed three-term recurrence relations in Appendix B. 

In contrast to scalar field case, the even and odd sequences are coupled 
each other through the non- vanishing S n terms. Dividing Eq. (|20|) for n = 3 by 



Eq.(20) for n = 5, we get the relation between the ratio of the first two terms 



of odd sequence and that of even sequence as follows, 

03 = "5 jf +75 + £5 Q4 
«2§f + 73 + £30«2 

■^ 2 ' + ff7 2 ')+75 + e 57 ]f 



_ 104/271 fa 1 ^ tt4 q 4 

a 37i f+73-e3^ a 2 

where the second equality comes from Eq.(|l8|) and the transformed three-term 
relation ( |Sl| ) for n = 2. The quasinormal mode boundary conditions are sat- 
isfied when the continued fractions ( p3| ) and ( |34| ) satisfy Eq.(|3^). This is the 
characteristic equation for the quasinormal frequencies of the extremal Reissner- 
Nordstrom black hole. 



Eq.(35) can be shifted an arbitrary number of times, n, to yield an equation 
for the ratio of successive terms of higher n: 



Qn+6 , , a n +2 

dn+3 _ a »+ 5 a„ +4 T ln+5 -t- tn+5 q k+4 a » +4 

«n+l " «™+2g±f + 7n+3 + ^+3^ ' 



(36) 



Eq.(35) and Eq.( p6[) are completely equivalent since every solution to (|35|) is 
also a solution to ( |36| ) when we substitute a continued fraction into each ratio 
of a;/ai_2, and vice versa. Eq.(|36|) can be used as a check of results obtained 
from Eq.®. 
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3 Numerical Results 



The eigenvalue equations (g6|) , ( p7j ) and (|35|) derived in the previous chapter are 
solved using the HYBRD subroutine distributed in MINPACK libraries. The 
three least-damped modes of (a, I) = (2,2), (2,3), (2,4), (1,1), (1,2), (1,3), 
(0, 0), (0, 1), (0, 2) for the extremal Reissner-Nordstrom black hole are listed in 
Table 1. We compare our results with Leaver's results|ll| of nearly maximally 
charged case (99.99% charged), which should be divided by two because of the 
difference of our scaling from his. These modes appear as complex conjugate 
pairs in p because the eigenvalue equations are characterized by only real coeffi- 
cients. We also compare them with the third order WKB quasinormal frequen- 
cies of the extremal Reissner-Nordstrom black hole. The WKB quasinormal 
frequencies are obtained using the same formula of Kokkotas and Schutz]l4]|, 
who calculated the frequencies of charged black holes but did not show those of 
maximally charged black holes in the paper. 

Our results are in agreement with other results obtained by Leaver and those 
obtained through the WKB method within the accuracy of a few percent. The 
difference between WKB frequencies and ours is ascribed to the fact that the 
WKB method gives only the approximate values. Indeed, the tendency that the 
discrepancy grows with the mode number suggests the breakdown of the WKB 
approximation for higher n. The difference between Leaver's results and ours is 
from the numerical error caused by the breakdown of his series expansion at the 
limit of maximal charge. Thus we are sure that our quasinormal frequencies of 
the extremal Reissner-Nordstrom black hole are the most accurate. Our results 
and his are plotted in Fig.l. 

Notice that the quasinormal frequencies of (s, I) = (2, 2) and those of (1, 1) 
approach each other with increasing Q and coincide in the extremal limit. We 
show this curious coincidence in the limit of maximal charge explicitly in the 
Fig. 2. The trajectories of the first four lowest modes of the third order WKB 
gravitational quasinormal frequencies of I — 2,3,4,5 are plotted. These quasi- 
normal frequencies coincide with the frequencies of electromagnetic wave of 
I = 1, 2, 3, 4 in the limit of maximal charge, respectively. The first order WKB 



quasinormal frequencies are given by ui^ = V pea k — i{n + l/2)w — 2V„" a/£ using 



the Zerilli-Moncrief Potential V, where a subscript n is a mode number. The 
differences of V pea k and V" ak between V\ for I — 1 and V% for I = 2 are calcu- 



lated to AV peak /V peak ~ 2.79565 x 1CT 16 and AV'' eak /V'' eak ~ 5.46819 x KT 15 , 



respectively. This explains the coincidence of a few least-damped quasinor- 
mal frequencies of these two modes. Even for higher overtone modes n ~ 30, 
however, we numerically find a remarkable coincidence in the quasinormal fre- 
quencies of gravitational waves of I and those of electromagnetic waves of Z — 1 
for a wide range of multi-pole indices, 2 < I < 10, through our method. 

These two potentials, V\ with multi-pole index I and V2 with multi-pole 
index l + l, indeed give the same transmission and reflection amplitude for each 
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a, whose evidence is presented in Appendix C. It implies that the quasinormal 
frequencies for both potentials are identical. 

4 Conclusions and Discussion 

We improve the continued fraction method of computing quasinormal frequen- 
cies of charged black holes in order to determine the quasinormal frequencies 
in the extremal limit and solve our formulae of eigenvalue equations by using a 
root search program. Our results of quasinormal frequencies for the maximally 
charged black holes are consistent with any other previous results fl|, @ and 
are the most accurate of all. 

Our procedure is summarized as follows. We expand the solution about the 
ordinary point and get the recurrence relation of the expansion coefficients. In 
our method we have to examine the boundary conditions at both boundaries. 
Then, we need to divide the sequence into an odd sequence and an even sequence 
to examine convergence of both series. The convergence of the series at both 
sides occurs if and only if each sequence is minimal. The ratio of successive terms 
of each sequence is given by a continued fraction. Consequently, the eigenvalue 
equation is constructed by substituting the continued fraction into the ratio of 
successive recurrence relations. 

The applicability of our method is free from types of singularities of the 
equation, in contrast to Leaver's method which is not useful when the equation 
has two irregular singular points at both boundaries. Though we can use our 
method to determine quasinormal frequencies of other black holes such as the 
Schwarzschild, non-extremal Reissner-Nordstrom and Kerr black holes, we actu- 
ally did not follow this procedure in these black holes. For these black holes, his 
method is easier to obtain the quasinormal frequencies because the recurrence 
relation and the eigenvalue equation are simpler than ours. In our method, it is 
generally difficult to give the explicit forms of the recurrence relations of both 
even and odd sequences. In that case we will have to eliminate the even or odd 
sequence by numerical elimination steps. For the extremal Reissner-Nordstrom 
black hole, the recurrence relation, Eq. ( |2p| ) , is not so complicated in that it has 
no /3 n -term and furthermore <5 n -term has no dependence on n, which makes our 
procedure rather easy to accomplish. 

We believe our method can be also used for the extremal Kerr black hole 
which has an irregular singular point at the horizon in the equation for radial 
part. Some authors |l6| [l7| , [TH obtain the series of infinite numbers of quasi- 
normal frequencies of the extremal Kerr black hole, which accumulate onto the 
critical frequency. Leaver also finds a tendency of an accumulation of quasinor- 
mal frequencies when the Kerr parameter approaches the maximum value. It 
is interesting to investigate what will happen to the distribution of quasinor- 
mal frequencies in the limit of maximal rotation, where Leaver's method is not 
valid. We are in preparation for this problem through our method of computing 
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quasinormal frequencies. 

We accidentally find a curious coincidence in the quasinormal frequencies of 
gravitational perturbations with multi-pole index I and those of electromagnetic 
perturbations with I — 1. These two modes are completely decoupled but they 
have the same amplitudes of transmission and reflection. In that case only 
the difference is a phase shift of transmitted or reflected wave, that means the 
quasinormal frequencies for both modes are identical. The situation is very 
similar to the coincidence of the transmission and reflection amplitudes of odd 
parity perturbations with those of even parity perturbation |fl5| . At present 
there is no easy way to understand the hitherto unobserved coincidence but it is 
interesting that it occurs only in the extremal limit, where the black hole may 
have an unknown symmetry. 
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Appendix A 

Here we summarize a procedure to obtain a recurrence relation for the odd se- 
quence and that for the even sequence of a n . First we write down the recurrence 
relation, Eqs.(^8|), ( |l9| ) and (p0|), in the matrix expression and separate the odd 
sequence from the even sequence as follows, 
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where the term S n is constant therefore we write it as 5 here, which makes us 
able to accomplish the below procedure explicitly. If we substitute the right 
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hand side of Eq.(p8[) into the left hand side of Eq.(|: 
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the even sequence is eliminated and then we can obtain a recurrence relation 
including only b n which was presented in Eq. (^0|) and the recurrence coefficients 
arc explicitly given by 
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Similarly we get a recurrence relation for c„ and the recurrence coefficients 
are given by 



a n = 1 — 5n + An , 

(3 n = (-1 + 2n)(2 + 2 A - 8n - AAn + I2n 2 - an" 

-q s + 2nq s + 6p + A8np - A8n 2 p + 32p 2 - 64n/j 2 ), 

72 = 41 + 2AA + AA 2 - 12q s - AAq s + q 2 + 3Q6p + 96Ap - 
+948p 2 + 128/lp 2 - 6Aq s p 2 + 1536p 3 + 1024p 4 , 

7„ = 38 + 24 A + 4^ 2 + 24n 4 - 12q s - AAq s 
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(45) 
(46) 
(47) 
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-6Aq s p 2 - 1536p 3 + 1024p 4 + 32n 3 (-3 + 7p) 
+2n 2 (77 + 8A- Aq s - 336p + 432p 2 ) 

+4n(-29 - 8,4 + Aq s + 189p + 2AAp - I2q s p - 432p 2 + 384p 3 ),(48) 
6 2 = -8 + AA + 6q s - q 2 + Yip - YlAp + 6q s p - 76p 2 

-16 Ap 2 + 8q s p 2 - 336p 3 - 256p 4 , (49) 
S n = -42 -6A- 16n 4 - (4 - q s f + 3q s + 32n 3 (3 - bp) 
+402/3 + A8Ap - 2Aq s p - 1160p 2 - 32Ap 2 + 16q s p 2 
+ 1344p 3 - 512p 4 + 4n 2 (-51 - 2A + q s + 180p - 144p 2 ) 
+4n(45 + 6A- 3q s - 2A7p - 8Ap + Aq s p + A32p 2 - 224p 3 ), (50) 
e n = 21 - 76n + 83n 2 - 32n 3 + 4n 4 - 152p + 332np - 192n 2 p 

+32n 3 p + 332p 2 - 384np 2 + 96n 2 p 2 - 256p 3 + 128np 3 + 64p 4 . (51) 



Appendix B 

We give the Gaussian elimination steps to transform the five-term recurrence 
relations into three-term recurrence relations. For b n the step defined by 
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for n > 3 , (59) 

transforms the five-term recurrence relation into a four-term recurrence relation 
and the second elimination step below, 
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brings us a final three-term recurrence relation in the following, 
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For c n the successive elimination steps in the following, 
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(67) 
(68) 
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(70) 
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for n > 3 , (74) 



(75) 
(76) 
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< = a' n , (78) 
K = K-<-iKK-i, (79) 

In = in-Pn-lS'n/in-1, &» 71 > 2 , (80) 

give us a three-term recurrence relation, 

<c„+i + BnC n + 7"c«-i - 0, for n > 2 . (81) 

Appendix C 

Now we summarize here the relation between Vi for multi-pole index I and V2 
for multi-pole index 2 + 1. If we introduce a new function /, 

r - 1 

/ = — > (8 2 ) 

which vanishes at both boundaries, r — > 1 and r — > 00, then two potential are 
given by 

Vi = + (/ + l)-g--4/ 3 + (Z + l) 2 / 2 , (83) 
F 2 = -(Z + l)-^--4/ 3 + (2 + l) 2 / 2 . (84) 

GET** 

If these two different potentials yield the same amplitudes of transmission and 
reflection, the following integrals for two potentials should be the same; 



h 



= j Vdr*, (85) 
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h = J V 2 dr«, (86) 

h = f(2V 3 + V ,2 )dr*, (87) 

h = J(5V 4 + 10VV ,2 + V" 2 )dr*, (88) 

h = J (14V" 5 + 7W 2 V 12 + UVV" 2 + V"' 2 )dr^ (89) 

h = 



Due to the simple forms, Eq.(|83|) and Eg. (|84|) , we can see the coincidence of these 
integrals, I\ ~ Jio, for two potentials. The explicit forms of above integrals, 
l\ ~ Is, are given by 

7i = (1 + 6Z + 3Z 2 )/3, (90) 
J 2 = (4+18Z + 93Z 2 + 84Z 3 + 21Z 4 )/630, (91) 
7 3 = (23 + 120/ + 476/ 2 + 1560Z 3 

+1820Z 4 + 858Z 5 + 143Z 6 ) /45045 (92) 
h = (1648 + 8736/ + 327162 2 + 105868Z 3 

+252567Z 4 + 316540/ 5 + 203490/ 6 

+64600Z 7 + 8075? 8 )/9399380, (93) 
h = (25160 + 125820? + 437634Z 2 + 1336860/ 3 

+3430245Z 4 + 6563970Z 5 + 8099427Z 6 + 6085800Z 7 
+2687895Z 8 + 642390Z 9 + 64239Z 10 )/1003917915. (94) 
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s = 2, 1 = 2 


s = 2, / = 3 


s = 2, 1 = 4 


11=0 

Leaver 

WKB 

n=l 

Leaver 

WKB 

n=2 

Leaver 

WKB 


(-0.083460,0.43134) 

(-0.083645,0.43098) 

(-0.08349,0.43013) 

(-0.25498,0.40452) 

(-0.257055,0.39309) 

(-0.25675,0.40076) 

(-0.44137,0.35340) 

(-0.442035,0.353515) 

(-0.44210,0.35136) 


(-0.085973,0.70430) 

(-0.08596,0.70398) 
(-0.25992,0.68804) 

(-0.26014,0.68701) 
(-0.44007,0.65624) 

(-0.43986,0.65575) 


(-0.087001,0.96576) 

(-0.08700,0.96563) 
(-0.26212,0.95381) 

(-0.26218,0.95339) 
(-0.44064,0.93020) 

(-0.44056,0.93004) 




s = 1,/ = 1 


s = 1,1 = 2 


s = 1,1 = 3 


n=0 

Leaver 

WKB 

n=l 

Leaver 

WKB 

n=2 

Leaver 

WKB 


(-0.083460,0.43134) 

(-0.08343,0.431415) 

(-0.08349,0.43013) 

(-0.25498,0.40452) 

(-0.259705,0.40602) 

(-0.25675,0.40076) 

(-0.44137,0.35340) 

(-0.44260,0.35347) 

(-0.44210,0.35136) 


(-0.085973,0.70430) 

(-0.086205,0.704075) 

(-0.08596,0.70398) 

(-0.25992,0.68804) 

(-0.26256,0.68315) 

(-0.26014,0.68701) 

(-0.44007,0.65624) 

(-0.4408,0.655675) 

(-0.43986,0.65575) 


(-0.087001,0.96576) 

(-0.08700,0.96563) 
(-0.26212,0.95381) 

(-0.26218,0.95339) 
(-0.44064,0.93020) 

(-0.44056,0.93004) 




s = 0,1 = 


8 = 0,1=1 


s = 0,1 = 2 


n=0 

WKB 

n=l 

WKB 

n=2 

WKB 


(-0.095844,0.13346) 

(-0.10371,0.12109) 

(-0.33065,0.092965) 

(-0.33742,0.09157) 

(-0.58833,0.075081) 

(-0.57164,0.05056) 


(-0.089384,0.37764) 

(-0.08936,0.37570) 

(-0.27614,0.34818) 

(-0.27828,0.34392) 

(-0.48643,0.29846) 

(-0.48145,0.29661) 


(-0.088748,0.62657) 

(-0.08873,0.62609) 

(-0.26909,0.60817) 

(-0.26944,0.60677) 

(-0.45820,0.57287) 

(-0.45750,0.57254) 



Table 1: The computed quasinormal frequencies are listed. Our results are very 
consistent with the values of WKB quasinormal frequencies. The relative errors 
of WKB complex and real gravitational quasinormal frequencies for I = 2 from 
our results are (0.04%, 0.03%), (0.07%, 0.1%) and (0.1%,0.5%) for n = 0,1,2, 
respectively. Leaver's value has rather large discrepancy of 3% for the real 
frequency of n = 1 gravitational wave. 
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Figure 1: Quasinormal frequencies for I = 2 gravitational wave and I = 1,2 
electromagnetic wave are plotted in the w-plane with the results of Leaver [13]. 
Each solid line is a trajectory of the first three lowest I = 2 gravitational mode, 
parameterized by the charge Q. It has a tendency to coincide at the right end 
with I = 1 electromagnetic mode which is shown as a small dashed line. Marks 
from left to right correspond to Q = 0, 0.4, 0.8, 0.9999 quasinormal frequencies 
of Leaver. The frequencies of extremal black holes we computed are plotted as 
diamonds. Dashed lines arc / = 2 electromagnetic modes. 
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Figure 2: Solid lines and dashed lines are trajectories of the third order WKB 
quasinormal frequencies of the gravitational and electromagnetic wave, respec- 
tively. Each left endpoint of lines corresponds to the quasinormal frequency 
of a charged black hole of Q — 0.8, and each right endpoint corresponds to 
the frequency in the limit of maximal charge. A trajectory of the gravitational 
quasinormal frequencies depicted in a solid line meets at the right end with a 
corresponding dashed line, which is a trajectory of the electromagnetic quasi- 
normal frequencies belonging to lower multi-pole index by one. 
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